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Abstract 

We study numerically the dynamics of quantized vortices and solitons 
induced by an oscillating potential inside a trapped Bose-Einstein con- 
densate. The dynamics of the topological defects is much different from 
the case for a linear uniform object; the metamorphosis between vortices 
and solitons is characteristic of the dynamics. We discuss how vortices 
are nucleated by an oscillating potential. 

PACS numbers: 03.75.Lm, 03.75.Lm, 03.75.Kk 

1 Introduction 

One of the most difficult problems in classical physics is turbulence. It is gener- 
ally considered that turbulence consists of vortices. The circulation of vortices 
has an arbitrary value, and the vortices repeat nucleation and annihilation. 
There is a similar phenomenon in a quantum fluid, which is called quantum 
turbulence (QT). In contrast to classical turbulence (CT), the circulation of vor- 
tices in QT takes a discrete value, and quantized vortices are more stable because 
they are topological defects [1]. Thus the constituent elements of QT are more 
definite than those of CT. The Kolmogorov law, which is the most important 
statistical law in turbulence, is confirmed in QT numerically with the Gross- 
Pitaevskii (GP) model[2]; these observations show some significant similarity 
between CT and QT in spite of the difference of the nature of vortices. 

A soliton is also important for quantum fluids, and there exists much liter- 
ature on the subject. Solitons and vortices are the representative topological 
defects in quantum fluids. Then it would be significant to study the synergy 
dynamics between vortices and solitons. 

Quantum turbulence has been studied in superfluid *IIe and '^He for these 
years [3]. Recently, QT is realized also in an atomic Bose-Einstein condensate 
(BEC) [4]. The key word for common technique making QT experimentally is 
oscillation. For example, QT in He is realized by oscillating the objects such as 
wires, grids, spheres and tuning folks, while oscillating a trapping potential leads 
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to QT in atomic BECs. However, there are not so many works with oscillating 
potential inside trapped BECs[5, 6]. 

In this paper, we study synergy dynamics of solitons and vortices induced 
by an oscillating potential inside a trapped BECs, and discuss the vortex nucle- 
ation characterized by oscillation. This approach should open a new aspect of 
quantum hydrodynamics. 

2 Formulation 

We consider a dilute atomic gas BEC with a macroscopic wave function, assum- 
ing that the condensate is pancake shaped. Then the macroscopic wave function 
tp{x,t) obeys the two-dimensional Gross-Pitaevskii equation: 

d 

i%—i){x,t) = - — A^{x,t) + (1) 

where m is a particle mass, g is an interaction parameter for two-dimensional 
case, and V{x,t) is potential. The macroscopic wave function ip{x,t) is nor- 
malized by the total particle number N . We suppose that the condensate is 
trapped by a harmonic potential and penetrated by a Gaussian potential, so 
V{x,t) = 14rap(a:) + K,bject(a;,i)- Here Vua.p{x) is harmonic as usual, T4rap(a;) = 
\m{ijS^x'^ -t-w^y^), and V^bject(!E,f) is represented by 

K,bject(r,t) = yoexp[-((ar - XQ{t)f + y'')/d% (2) 

where (t) is the x-coordinate of the center of the potential and d is the radius 
of the Gaussian potential. We move the Gaussian potential like xo{t) = €sm{ujt) 
with the oscillating amplitude e and the frequency u). The velocity of the po- 
tential is Vp = euj. 

We numerically solve Eq. (1) to investigate the dynamics of the topological 
defects (vortices and solitons). Neither dissipation nor noise is included in our 
calculations. 

3 Numerical results 

We choose the parameters as g = 4.19 x 10^''^J/m^, m = 1.42 x 10^^ kg, 
= 6.6 X 10"', ujx = 2tt X 5 /s, ujy =2-!t x 25 /s, d = 6fj,m, e = 7/im, u = 100 /s 
and Vo = GO^no, where no is the density near the center of the condensate. In 
these parameters the Thomas-Fermi radii are Rx = 59.6 /im and Ry = 11.9 /im, 
and the sound speed near the center of the condensate is 1.33 x 10^/xm/s. We 
start the simulation with an initial stationary state of Fig. 1 obtained by the 
imaginary time step of the GP equation. 

This section consists of four parts. We describe elementary dynamics for 
vortices and solitons in the first three parts and show the global dynamics in 
the last part. One of characteristics of dynamics is metamorphosis between 
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Figure 1: (color online) Initial state of the condensate density: The x and y sizes 
of the box are 140 /im and 33 /im. Thomas- Fermi radii are — 59.6 /im and 
Ry = 11.9 /im. The low density region in the center comes from the Gaussian 
potential with the radius 6 /im. The right bar shows the range of the dimen- 
sionless density. This bar is common through Figs.l 6. 

vortices and solitons. The solitons appearing in our calculations are related 
with rarefaction pulse [7, 8] which is finite amplitude sound wave. 

Here we discuss heating in this system because the oscillating potential may 
increase the temperature and break a BEG. We can numerically calculate the 
increase of the total energy, estimating a change of temperature by using specific 
heat for equilibrium state. The increase of the temperature is found to be the 
order of 10 nK. Hence the heating can be ignored, and our numerical calculations 
with the GP model are valid. 

Nucleation of a vortex pair 

We discuss the nucleation of a vortex pair. Two kinds of vortices appear in 
our numerical calculations. One is a usual vortex which wears the profile of the 
condensate density, and the other is a ghost vortex [9] which exists in the very 
low density region and cannot be observed by the density measurement. We 
consider nucleation of both kinds of vortices in this work. Then there are two 
critical velocities; Vci for a pair of ghost vortices and Vc2 for a pair of vortices 
with Vci < Vc2- When the velocity v-p of the potential is smaller than Vci, 
any vortices are not nucleated in the condensate. As v-p is increased to exceed 
Del, a pair of ghost vortices is nucleated inside the Gaussian potential. If Vp 
is smaller than Vc2, however, ghost vortices are immediately annihilated. Thus 
ghost vortices repeat nucleation and annihilation inside the potential, so that 
they cannot get away from it and vortices are never nucleated in the condensate. 
When Vp exceeds Vc2^ the vortices are nucleated eventually. 

This process is shown in Fig. 2. Firstly, the velocity field like back-flow 
is induced by the potential in Fig. 2(a). Secondly, a pair of ghost vortices is 
nucleated inside the potential in Fig. 2(b). Thirdly, the ghost vortices get away 
from it, then a vortex pair is nucleated in Fig. 2(c). The repeated nucleation and 
annihilation of a pair of ghost vortices are characteristic of the oscillating case 
with Vci < Vp < Vc2- In the present case, Vd = 500/im/s and Vc2 — 580 //m/s, 
but Wci and Vc2 strongly depend on the amplitude e. The critical velocity Vc2 
is smaller than the sound speed 1.33 x lO'^/^m/s, similar things are reported by 
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(a) (b) (c) 

Figure 2: (color online) Vortex nucleation: (a) t=0.00253s, (b) i=0.00490s, (c) 
t=0. 0118 s. The upper and lower pictures denote the density and phase profiles. 
The X and y sizes of box are 36 fim and 33 /im through Figs .2 5. The symbol 
+ denotes a vortex with the clockwise circulation and - denotes a vortex with 
the counter-clockwise circulation. Arrows refer to the direction of velocity of 
the potential. 

other authors. The dependence should be systematically investigated, reported 
soon elsewhere. 

Dynamics of vortices 

In this subsection, we show the dynamics of vortices in Fig. 3, immediately after 
the potential starts to oscillate. The dynamics of vortices by an oscillating po- 
tential is quite different from that by a linear uniform moving potential. Figure 
3 shows the characteristic dynamics of vortices. As the potential moves, a vor- 
tex pair is nucleated behind the potential through the process in Fig. 2, and has 
the impulse whose direction is the same as the velocity of the potential in Fig. 
3(a). Consequently, the pair follows the potential. After the oscillating poten- 
tial changes the direction in Fig. 3(b), another pair of ghost vortices is nucleated 
inside the potential and the old pair collides with the potential. Through this 
collision, the ghost vortices come out of the potential to become a new vortex 
pair, which reconnects with the old pair in Fig. 3(c), then the rearranged pairs 
leave the potential and move towards the surface of the condensate in Fig. 3(d). 
This phenomenon is peculiar to the oscillating potential (AC), not observed in 
linear moving potential (DC). Note that the ghost vortices are nucleated in the 
surface of the condensate as the potential moves. These pairs in Fig. 3(d) arrive 
at the surface, and reconnect with ghost vortices. Afterward, new vortex pairs 
move along the surface of the condensate and head toward the bow of it [10]. 
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Figure 3: (color online) Dynamics of vortices: (a) i=0. 0507 s, (b) t=0. 0551s, 
(c) t=0. 0571s, (d) t=0. 0596 s. These pictures denotes the density profile. The 
symbol + denotes a vortex with the clockwise circulation and - denotes a vortex 
with the counter-clockwise circulation. Arrows refer to the direction of velocity 
of the potential. 
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Figure 4: (color online) Nucleation of solitons: (a) t=0. 0929 s, (b) t=0.106s, (c) 
t=0.111s, (d) t=0.114s. These pictures show the density profile. The symbol + 
denotes a vortex with the clockwise circulation and - denotes a vortex with the 
counter-clockwise circulation. Arrows refer to the direction of velocity of the 
potential. In the process from (a) to (b), reconnection of vortices occurs, then 
two new vortex pairs, which will become solitons, are created. These pairs are 
depicted by closed loops with broken lines. The closed loops with solid lines in 
(d) show the solitons. 

Dynamics of solitons 

Our calculation shows the interesting and remarkably dynamics that solitons 
are nucleated by annihilation of vortex pairs and the collapse of solitons makes 
vortex pairs. 

The annihilation of vortex pairs makes solitons. Figure 4 shows this process. 
The vortex configuration in Fig. 4(a) is obtained by repeating the processes in 
Fig. 3 two times. Then the vortex pairs rearrange near the surface in Fig. 4(b) 
and new vortex pairs move to the central part of the condensate in Fig. 4(c). 
These new vortex pairs are depicted by closed loops with broken lines. Then 
the annihilation of the vortex pairs nucleates solitons in Fig. 4(d). The soliton 
created by the annihilation are related to rarefaction pulse. 

Next we show that the collapse of the solitons leads to the nucleation of 
another vortex pair. This process occurs after Fig. 4, shown in Fig. 5. After Fig. 
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(a) (b) (c) (d) 

Figure 5: (color online) Collapse of solitons to vortex pair and their dynamics: 
(a) t=0.126s, (b) <=0.131s, (c) t=0.136s, (d) i=0.138s. These pictures show the 
density profile. The symbol + denotes a vortex with the clockwise circulation 
and - denotes a vortex with the counter-clockwise circulation. Arrows refer to 
the direction of velocity of the potential. The closed loops with solid lines in (a) 
show the solitons. The symbols ® and denote the vortices created through 
the collapse of the soliton. 




Figure 6: (color online) Trajectories of vortices: The x and y sizes of the box 
are 140 /xm and 33 /im. The snapshot is the density profile at t=0. 0845 s. The 
symbol + denotes a vortex with the clockwise circulation and - denotes a vortex 
with the counter-clockwise circulation. Arrows refer to the direction of velocity 
of the potential. The curves show the trajectories of vortices which do not change 
to the solitons. The same motion occurs in the left side of the condensate. 

4(d), the solitons go through each other without changing the shape, which is 
characteristic of soliton. Then Fig. 5(a) is obtained. The solitons go toward the 
surface in Fig. 5(a) and the solitons collapse to vortex pairs [11] near the surface 
in Fig. 5(b). The reason why this collapse occurs in the surface may be that the 
phase of the macroscopic wave function is disturbed in the low density region. 
Later, the new vortex pairs reconnect with the ghost vortices in Fig. 5(c) and 
the new pairs head toward the bow of the condensate in Fig. 5(d). 

Global dynamics of vortices 

Global dynamics of vortices consists of the above elementary processes. In Fig. 
6 we confine our description to right region of the condensate. The vortices 
nucleated by the Gaussian potential leave it by the process shown in Fig. 3. 
Subsequently, the vortices reach the surface and reconnect with ghost vortices. 
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The resulting pairs at the upper and lower sides of the condensate go to the 
bow along the surface, then the pairs meet near the bow, and again reconnect 
with each other. As a result appear a pair of vortices which heads toward the 
center of the condensate. 

4 Conclusion 

We perform the numerical research with the GP model to investigate the dynam- 
ics and nucleation of vortices and solitons induced by an oscillating potential in 
a trapped BEG. Our calculations found three interesting phenomena. Firstly, 
the dynamics of AG case is much different from that of DG case in many re- 
spects, while they are simmilar in some respects. For example, the trajectories 
of the vortices shown in Fig. 6 are observed in a DG experiment [12], but the dy- 
namics shown in Fig. 3 is characteristic of the AG case. Secondly, the dynamics 
and nucleation of vortices are deeply involved with the ghost vortices. In par- 
ticular, we consider that the ghost vortices may be essential to the nucleation 
for the vortices. Thirdly, the synergy dynamics between vortices and soliton is 
realized. This is related with the rarefaction pulse. One can enjoy the movie of 
the numerical simulation in our web site[13]. 

The parameters used in our calculations are realistic in experiments of atomic 
BEGs. Hence, the dynamics obtained by our study can be observed actually in 
experiments. We consider that oscillating Gaussian potential may lead to QT 
so that this method may be a new technique for making QT in atomic BEGs. 
The details reviewed will be reported soon elsewhere. 
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